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Abstract. We derive simple linear, inliomogeneous recurrences for the variance of the index by utilising the 
fact that the generating function for the distribution of the number of positive eigenvalues of a Gaussian unitary 
ensemble is a T-function of the fourth Painleve equation. From this we deduce a simple summation formula, 
several integral representations and finally an exact hypergeometric function evaluation for the variance. 



1. Introduction 

The exact characterisation of distribution functions in random matrix theory using integrable systems 
theory is very prevalent throughout random matrix theory As well as bringing to the fore rich mathematical 
structures, such characterisation has the practical consequence of both providing efficient computational 
schemes for the distributions, and in the determination of their asymptotic properties. 

Here, inspired by a recent work of Majumdar et al |[T4|, we will undertake a study from the viewpoint of 
integrable systems of the so-called index distribution for the n X n Gaussian unitary ensemble (GUE). The 
index distribution is a discrete probability distribution for the number of positive eigenvalues. It is therefore 
specified by the probability p(n+, n) for there being n+ positive eigenvalues out of the n eigenvalues in total. 
We refer to the Introduction of IIT4ll for an extended discussion as to the interest in the index distribution for 
its relevance to studies in disordered systems, landscape based string theory and even quantum cosmology. 

Explicitly, our interest is in the integrable systems context of the variance of the index, 

n 

(1) A„ := {{n^ - {n^)f) =Yj{n^- {n^)f p{n^, n). 

n+=0 

Note that by the symmetry M+ i-^ n -«+ we must have the average number of positive eigenvalues, (m+) = 
Our working expression, see Eq. (137) of [MJ, is obtained by introducing the "Poissonisation" of the index 
distribution 

n 

(2) Zn{s):=J^e-'"yin^,n). 

n+=0 

Most importantly, making use of the explicit form of the joint eigenvalue probability density function for 
the GUE, this has the multiple-integral representation 



(3) Z„(s) = — dAi... dA„ e-^;^>''H^M J] (A^ -A,c)^ n > 1, Zq = 1, 

where 



l<j<k<n 



[1, A>0 

(4) 0(A) = U, A = 0. 

[o, A < 

is the Heaviside step function. The normalisation N„ is the Gaussian weight form of the Selberg integral 
(see Proposition 4.7.1 of [TO]) 

dAi... dAne-^''"^ Yl {A,~Akf, 

l<j<k<n 



and has the evaluation 

(6) N„ = 2-2"("-i'7t"''2m!G(« + 1), 

where G(z) is the Barnes G function (see of 5.17 of [22]). From the generating polynomial ^ we can reclaim 
the probability itself and this is given by 

(7) p(«+, n) = ^1^" J '^'^^ ■ • ■ X '^^"^ L '^^"'^^ ■ ■ ■ X ^' n ^'^^ " ^ ^ 

Using the form l|2ll we seek to compute the variance according to the formula 



1 

(8) A„ = ^JjZn 



- ln\ n>0. 



s=0 

From the integrable systems viewpoint, the significance of l|2]l is that it can be written in the Hankel 
determinant form (see HOJ ), Eq. (5.75)) 



(9) 



Z„(s) = ^ det[ dA A'^''e-'''-''^'%^j,,^„.^, n>l. 



The work fT4] used ^ to derive various summation and an integral formula, and to deduce the leading 
order asymptotics of A„. However no tie in with integrable systems theory was made. In fact Z„(s) is a 
T-function of the fourth Fainleve equation and much of the fundamental theory for this system has already 
been worked out in a study by the present authors [llj. In particular simple and effective formulae to 
characterise the variance of the index will be derived as an application of the results in this latter study. 
We also remark that the techniques employed here can be generalised to treat the higher moments of the 
index distribution with a little more labour, and in fact can be applied to other jS = 2 ensembles where the 
T-function theory has been similarly developed. 

2. The GUE and the fourth Fainleve Equation 

2.1. Okamoto x-f unction theory. We recount some of the fundamental facts in Okamoto's T-function theory 
of the fourth Fainleve equation |2T|| , |20| , as given in the reformulation by Noumi and Yamada flTl (see also 
the monograph lllSil ). As mentioned in the Introduction this theory was summarised in IfTTl and applied to 
the Gaussian unitary ensembles, and it is the key results given there that will be of the greatest utility in 
addressing our problem. 

The first key result is the symmetric form of the fourth Fainleve equation | 



Proposition 1. The fourth Fainleve equation is equivalent to the coupled set of autonomous differential equations 
(where ' = d/dt) 

/o =/o(/i-/2) + 2ao, 

(10) /; =/i(/2-/o) + 2ai, 

/2 = /2(/o-/i) + 2a2, 

subject to the constraint 

(11) fo+fl+fl = lt, 

with the fourth Fainleve transcendant given by y = -fi and where the parameters aj e IR with ao + ai + a2 = 1 are 
related to the conventional parameters 

(12) a = a^- a2, p = -2aj . 

The parameters {aQ,ai,a2) form a triangular co-ordinate system system in the root system attached to 
the extended affine Weyl group A^^^ and the fimdamental Weyl chamber given by < ao, fli, «2 < 1- 

The second feature of the fourth Fainleve system is the existence of a Hamiltonian formulation, which is 
essential to our theory. 
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Proposition 2. 1211 [131 The PIV dynamical system is a non-autonomous Hamiltonian system {q, p; t, H] with the 
Hamiltonian 



(13) H={2p-q- lt)pq - la^p - a2q, 

(14) = 5/0/1/2 + a2 /i - ai fi, 
and canonical variables q,p 

(15) q = -h, P = \f2- 

This allows us to define the x-function T{t) in terms of the Hamiltonian H{t) by 

(16) H(0=:|logT(0. 

We will construct a countable sequence of classical solutions to the fourth Painleve equation with parameters 
{uq + n,a\,a2 — n), indexed by n e Z>o. The seed solution, n - 0, locates the parameters on one of the walls 
of the affine Weyl chamber, a2 ~ 0, and this is one of the necessary conditions for a classical solution of 
the fourth Painleve equation. Thus we construct a sequence of T-functions t„(0 (which were denoted by 
T3[n]{t; a\) in Ull) that also satisfy the Toda lattice equation 

(17) ;glog(.-'>-«)T„) = C:^^. 

Proposition 7 of flU states, in part, that such x-fimctions are also Hankel determinants. 

Proposition 3. A sequence of solutions to the Toda lattice equation lfT7l > with parameters {uq + n,ai, -n), n > 0, 
starting from the Weyl chamber wall a2=Q is given by the determinant form 

(18) T„{t;ai) = Cndei\ ( {t - xy'^'^k-""^ dx] , n>l, To = 1, 

L U,j=Q,...,n-l 

where C„ is a constant independent oft and a\. 

However this result can be easily generalised in a way that is of direct significance to our problem. 

Proposition 4. A sequence of more general solutions to the Toda lattice equation W7\ with parameters (o-q + n, a\, -n), 
n > 0, starting from the Weyl chamber wall a2=0 is given by the determinant form 



(19) 



r„(t;ai;.5) = C„det[|j^ +£, J (f - x)-"i+'+^e-'^' n > 1, To = 1, 



where ^ is an arbitrary complex parameter and C„ is also independent ofE.. 

Proof. We do not enter into the details of the proof of this as the proof of Proposition 7 in fTTll carries over 
in this case in its entirety. □ 

Finally we construct pure recurrence relations in n for our dynamical system which can be shown to be 
equivalent to the second order difference equations of the first discrete Painleve equation dPi (see Proposition 

iiof im). 

Proposition 5. The symmetric variables {fj,n{t', «i)}/=o,i,2 satisfy a coupled system of first order recurrence relations 
inn>0 

In 

(20) /o,„ + /o,„-i =lt- /2,„ + — , M > 1, 

J2,n 

(21) /2,„+i+/2,„=2t-/o,„ + ^^^y^, n>0, 

/0,Ji 

together with fo^„ + /i,„ + /2,„ = 2t. 



In IITTll the symmetric variables were denoted by !/;[«]};=o,i,2- A number of relations enable one to recover 
the symmetric variables from the T-functions and vice-versa. Thus one has 

(22) f,„f,„=-2n + ^^^^, n>l, 

T 

which motivates the definition 

(23) gni^) := fiAmA^y 

The symmetric variables can be recovered using 
.24^ r2 ^ (g„ +2n)g„+i +2wg„ 

gn+2n{gn + 2n)g„.i+2ngn' 

(25) ^2 ( 

^ ' ' {gn+2n)g„^i+2ng„' 

modulo the ambiguity associated with the choice of the branch of the square-root. 
Also given in the work fTTll is a recurrence relation for the x-function directly. 

Proposition 6. A single recurrence relation for the T-functions is given in | |TT| hy Eq. (2.82), which states n>2 

(26) [t„+it„-i - 2«T^] [t„+it„-i + 2(ai - m)t^] 

X [t„+2T„t„-i - 2Tf,^^T„-i + 4M(ai - n)T„+iT^] [t„-2T„t„+i + 2t^_jT„+i + 4«(ai - m)t„-iT^] 

= ^lln+l'ln-l + T„+2T^_i [t„+iT„_i + (2ai + 2 - 4n)T^] + [t„+it„_i + (2ai - 2 - 4n)T?,] 

- \6n^{ai - nfjl + 16n{ai - n){ai - 2«)Tf,T„+iT„_i - 4(2n^ - 2ain + 1)t„t^^^t^_i| . 

Proof. In ITTll the /i, /2 system of coupled recurrences can be eliminated in favour of the T-functions through 
the relation 

(27) 4f2c„(b„-i - Cn)(b„ - c„) = {b„b„-i - clf, 
where the auxiliary variables are 

(28) bn ■= fl,nf2,n + fl.n+lfl.n+l + 2ai, 

and 

(29) c„ :=/i,,,/2,„—- — —— . 

f\,n J2,n + 2n 

Eq. l|26l l results from the substitution of (|22] | into this system. □ 

Our first key result is the following identification. 

Proposition 7. The moment generating function ^ is the n-th member of the classical x-function sequence of the 
fourth Painleve equation with independent variable t = 0, parameters (uq + n, a\, -n) with ao = 1, ai = and initial 
conditions £. = e~^ by 

(30) Z„{s) = ^^^^[^_^j „{Q;Q;e-'), n>0, 

where N„ has been given by With these specialisations the z-functions satisfy a simplified recurrence inn>2 

(31) T„+2Tf,T„-2 + T„+2Tf,_i [t„+iT„-1 - (4« - 2)t?,] + T„-2Tf,^i [t„+iT„-1 - (4« -h 2)t^] 

-A{2n^ + l) t2^iT„t^_i + 32«\„+iTf,T„-i - 16«*Tf, = 0. 

or alternatively the g-variables satisfy the recurrence 

(32) gl = {2n + gnfign + gn-l){gn + gn+l), n > 1. 
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Proof. This follows from a comparison of ^ and l fT9ll and the only issue that remains is to determine the 
pref actor C„. This can be done for 4 = 1 as C,, only depends on n and not t, f etc. Thus we have 



(33) 



T„(l) 



Cn2-i"^"-^^n"'^G{n + 1), 



where G(z) is the Barnes G function (see of 5.17 of [22]), which we substitute into (|3TJ. Noting that this only 
involves the second-ratio of C„ we make the ansatz 



(34) 



Cn+2 = C 



Cn ' 



and compute that the recurrence contains an overall factor of (c - 4)^. Thus we set c - 4 and take the simple 
solution C„ = 2"("-i). □ 

Thus our key quantity of interest is now expressed by the x-fimctions through 



(35) 



A„:= 



T„(^) 



n n 

- -r + TT/ « ^ 0- 
4 2 



Remark 1. Clearly t„{£,) is a polynomial in £, with degree degc ~ n, however the independent variable in our 
context is not a "time" variable in the conventional dynamics of the fourth Painleve system and does not 
appear in any of the characterising differential or difference equations. This parameter arises as a boundary 
or initial condition and will also be part of the monodromy data of the associated isomonodromic system. 
Our system is a classical one and therefore the monodromy is expected to trivialise in some way, i.e. the 
monodromy matrices will be upper or lower triangular or ±1 multiples of the identity, yet it will be a 
one-parameter system through ^. 

Remark 2. The seed solution is located at the corner of the Weyl chamber cto = ^,ct\ ~ 0, a2 ~ 0, and this is 
known as one of the conditions for the rational solution that generate the family of "generalised Hermite 
polynomials" IflTl , lU, [71, and are expressible as Schur functions for rectangular Young diagrams. However 
it is not clear there is any link between those polynomials and the ones treated here. 

Remark 3. In our identification. Prop. [71 we set ai = 0, however all of the analysis we will present carries 
through with ai non-zero. The case where such a situation is of relevance to our original problem is when 
«! is a negative integer and we are conditioning a fixed number of eigenvalues to lie at the origin. 

2.2. Low n Solutions. We present some data for the low n examples of our system, in order to furnish initial 
values in the n recurrences we are about to employ and also to furnish independent checks on our results. 

With the initial values for the x-functions T-i{£,) := 0, To(<^) := 1 and their Hankel determinant form for 
n>l 



(36) 



(37) 



T„(£) = 2"(«-i) det 



dx{-xy^'e 





- 0<i,j<n-l 



2« det [(5 + (-l)'^Or(i(f + ; + 1))]^^, 



y<"-l 



we easily compute the first few cases of the T-fimctions 



Ti(4) = - V^(5 + 1), 
T2(5) = ^ [71(5 + 1)2 -2(^-1)2], 
T3(5) = 2 V^(4 + 1) [71(5 + 1)2 - 3(£ - 1)2] , 

T4(£) = 8 [en^ii + If - 29n{e - if + 32(£ - 1)^] , 
T5(5) = 128 Vn(4 + 1) [727i2(£ + If - 4357i(£2 - \f + 656(£ - 1)*] , 
T6(5) = 4096 [4320ti3(4 + 1)^ - 337237x2(4 - lf{^ + 1)^ + 841687t(<5 - lf{E + 1)2 - 65536(4 - 1)^] ■ 



Correspondingly we compute the initial auxiliary variable g,, as 

goiQ = 
4(5-1)2 



gm = 



g2{i) = m - If 



[71(4 + 1)2-4(5-1)2] 
[7l(£ + l)2-2(5-l)2]2' 



(5 - 1)2 l^nHi + 1)" - 36ti(.£2 _ ^^2 + 64(4 - 1)4] 



n{£. + 1)2 [7^(4 + 1)2 _ _ 1)2]2 

The initial values of the symmetric variable /i „ and the leading terms of their Laurent expansions are 

2(4-1) 



(38) /i,o(4) 



V^(4 + l)' 



(39) =izi_(i:_lZ,(tzH!,oK-i)' 

Vtt 2 \n 4 y7T 
4(4 - 1)3 

(40) /i,i(4)= ^ 



V7l(4 + 1) [71(4 + 1)2 -2(4 -1)2]' 



(42) /„(«= 



V7l(4 + 1) [7t(4 + 1)2 - 3(4 - 1)2] [71(4 + 1)2 - 2(4 - 1)2]' 



(43) =i^-<izi£.(i^M^,0(4-l)« 

2^/n 4 ^n^' 

Those of /2,„ are 

(44) /2,o(4) = 0, 

[71(4 + 1)2-2(4-1)2] 

(46) =-pl^^^izl^o{^-lf, 

4 - 1 V" 

4V7i[4'-l][7z(4 + l)2-3(4-l)2] 

(47) /2,2(4)= ^ ^ 



[71(4 + 1)2 - 4(4 - 1)2] [71(4 + 1)2 - 2(4 - 1)2]' 

(48) =?<i_l)-<i_i)!,<!Ll|)|_i)!.o(i-i,., 

whilst the corresponding data of /o,„ are 

(49) 2(5-1) 



V^(4 + l)' 



(50) =-£_i,(iiiS-(i:iH!.o(4-i)«, 



_ V^(5 + l)[n(5 + l)2-4(5-l)2] 

^^^^ " - 1) [7l(.£ + 1)2 - 2(5 - 1)2] ' 

(52) =p^^^^izl^o{i-lf, 

(5 - 1) [7T(£ + 1)2 - 2(5 - 1)2] [571(5 + 1)2 - 16(5 - 1)^1 

(53) /o2(5) = — ^ 

V^(5 + 1) [71(5 + 1)2 - 4(5 - 1)2] [71(5 + 1)2 - 3(5 - 1)2] 

2 Vn 4 Vti 87T'^'2 

We remarked earlier that there is a sign ambiguity in deducing the symmetric variables from the T-functions. 
In computing the above variables we have adopted a particular and arbitrary choice for the sign of an initial 
variable, say that of /i^o as being +1. The recurrence system 1120121b fixes the sign of the other components 
at n = and all components subsequently for n > 0. 

Finally we give some initial values of the variance of the index itself 



Ai = 



1 

4' 
1 



A2 = r - 
2 n 



A3 = 7 
A4 = 

As = 
A, = 



3 _ ^ 

4 27t' 

1271 

5 145 
4 487i' 
3 1249 
2 ~ 32071 ' 



Ai8 = 

Al9 = 



9 1198597830455957 



A' 



20 



2 913593230950407Z ' 
19 22773358778663183 
T ~ 164446781571072071' 
^ _ 183365193212828149 
^~ 1249795539940147271' 



3. Linear recurrences for derivatives 
We recall that the full non-linear recurrence relation for the r-function, « > 2, is 

(55) Tn+2Tn-2Tl + In+l^l^i [(2 " 4«)t2 + T„-iT„+i] + Tn-lll+i [("2 " ^n)%l + T„-iT„+i] 

- Ibn^xl + 32m3t„_it3t„+i - 4(i + 2n^) tI_^t„tI^^ = 0. 

where the dependence on 5 is completely implicit, through the initial conditions. Our strategy is to 
successively differentiate this relation with respect to 5 and set 5^1 afterwards, utilising the results of 
earlier solutions at each step. This does not necessarily guarantee a linear difference equation but in fact we 
will find that this is the case. 

Proposition 8. The ratio 

T (1) 

(56) Xn = n>0, 

Th(1) 
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satisfies the linear recurrence relation 

(57) nXn^i - (n + 1)X„ - (n - 2)X„_i + (n - 1)X„_2 = 0, « > 2, 
which has the solution 

(58) X„ = |, n>0, 

g-fz^en f/ze znz'h'aZ conditions Xq = 0,Xi = j,X2 = 1. This implies the assertion given in the introduction that 



2 

(59) <n^) = -^T, 

T„ dt 



2 
4=1 



- -n, n> 0. 



Proof. We first set 

(60) T„(l) = 2"^"-^^'\"'^G(n + 1), n > 0, 
and define 

(61) T;,(l) = 2"("-W27^"/2G(n + l)X„, «>0. 

Differentiating l l55l l once and setting both t„(1) and t„(1) as given above leads to the result that the left- 
hand side is identically zero. Differentiating again and making the same substitutions gives the result that 
the coefficients of the second derivatives vanish identically however the first derivatives remain. After 
simplifying we find the relevant remnant factorises into a product of linear, difference equations 

(62) [(« - 1)X„_2 - (n - 2)X„_i - (« + 1)X„ + «X„+i] 

X [nXn-i -{n- 1)X„ - (« + 2)X„+i + (« + l)X„+2] =0, n>2. 

The second factor is the just the first under the replacement n ^ n + 1. A solution to the above difference 
equation could be one that satisfies the third order linear homogeneous difference equation l l57t for either 
the odd n or even n cases, or possibly both. The latter case applies here because the initial conditions enforce 
this. □ 

Continuing with this procedure we have the following result. 

Proposition 9. The variance of the index satisfies the third order inhomogeneous linear difference equation 

(63) nA„+i - (n + 1)A„ - (n - 2)A„_i + (« - 1)A„_2 = :^e„, n>2, 

with e„ as yet undetermined. 

Proof. When we differentiate (|55l l three times and make the previous substitutions for the T-function and its 
first derivative (including the solution given in | |58|| ) we find this identically vanishes. Let us set 

(64) t;;(1) = 2"'-"-^^'\"^^G{n + l)y„, n>0. 

Differentiating once more and making all of the above substitutions we can factorise the result into the 
product of two linear inhomogeneous difference equations 

(65) l-2n + 1 + 2(n - l)y„_2 - 2(n - 2)y„_i - 2{n + 1)Y„ + 2nY„^i] 

X [-2n - 1 + 2nYn-i - 2{n - l)y„ - 2{n + 2)y„+i + 2(n + l)y„+2] = 0, n>2. 

Again the second factor is just a copy of the first under n n + 1. However in this case there is an odd-even 
dichotomy and therefore only one of the factors vanishes - the other is an as yet undetermined n-dependent 
term which cannot be fixed by the recurrence relation (|55] l. It turns out that the even case vanishes and the 
odd is non-zero. Employing the change of variables 

(66) y„ = ^n{n - 2) + A„, « > 0, 

we have |[63] |. □ 

We can also carry out the successive differentiation of the g = fixfi recurrences (|32] | and deduce a system 
of linearised recurrences. 



Proposition 10. The first two Taylor coefficients of g„for alln>Oat^-l are given by 

(67) ^,,(1) = 0, 
and 

(68) g'„{l) = 0, 

whereas the second derivative is not uniquely determined by this method applied to (|32] I. 

Proof. In a pattern similar to the case of the T-function recurrences the odd derivatives vanish identically in 
a sequential way. The second derivative, for n > 0, gives 



(69) 



£-i(i)+^;,(i)][£(i)+^;,.i(i)]=o. 



The solution relevant here has g„(l) = ^„_]^(1) = for all n> Q. The fourth derivative gives a very similar 
result, namely 

(70) [£-i(i)+^:'(i)][£(i)+£.i(i)]=o. 

Subsequently we will see that g2,„+i(l) + Sim+i^^) = ^' however g'^Jl) + ^2,„+i(l) a 
Lemma 1. The inhomogeneous term is given by 

(71) ^'2m+l - kL(l) + 5W(1)' '"^o. 



Proof. If we successively differentiate the definition l|23t combined with the relation 
find, at the zeroth level 

(72) gn{l) = 0. 

At the next level we employ the definition (|56] l and its solution (|58] l to deduce that 

(73) ^:(1) = 0. 

This is consistent with our previous findings. At the second level we utilise (|64] l and 

(74) g;(l) = 2n[A„+i+A„_i-2A„], n>l. 
Then jTll follows immediately from (|63] |. 

The even and odd cases of A„ are related in our next result. 
Corollary 1. The odd and even cases of the variance of the index are related by 

2m + 1 ^ 

(75) A2„,+i = — A 



and set £, = 1 we 



to derive 



□ 



2m 



Im, 



m > 1. 



Proof. We observe that the difference equation | |63l l can be written as 
(76) 



1 - (-1)" 
rn - r„-2 = e„, n>2, 



where r„ := «A„+i - (n + 1)A„. In the even case of ((76t we have r2m = ''2m-2 = 
([75) follows. 



= ro = 0, and the relation 
□ 



Our recurrences have a structure allowing for a simple solution by discrete quadratures. 
Corollary 2. The variance of the index has the unique solution 



(77) 



i-(-i) 



^7(7 + 1) 2 



L 

1=1 



en+i 

71 



n > 0. 



Proof. The difference equation (|76l l can be solved by summation 

, w > 1, 



1 _/_iv' (^"^'^ 
(78) r„ = ^ 



eij+i + n 

/=1 

and we note the definition of r„ can be written in a telescoping form 

(79) = 

n + 1 n M(n + 1) 

and is also immediately soluble. In addition we utilise ri - -1/n. □ 
In order to evaluate the inhomogeneous term e„ we need to analyse the GUE fo, f2 coupled recurrences. 
Corollary 3. The inhomogeneous term e„ is given by 

1 r2(m + i) 

(80) e2,«+i--— 

TC^ T^(m + 1) 

Proof. We seek an expansion of GUE fo, f2 coupled recurrences separated into the even-odd sub-cases 

Am 

(81) foMQ + /0,2m-l(5) = -/2,2m(5) + 7 7TT, m > 1, 

j2,2in{C) 

(82) /0>,+l(«) + /0,2,„(4) = - /2,2m+l(«) + T^^^^, '« > 0, 

/2,2m+l(,Cj 

(83) /2,2m+l(5) + /2,2m(5) = " /o,2m(«) + ^ > 0, 

JQ,2m[Q) 

(84) /2,2m+2(«) + /2,2,„+l(«) = " /o,2,„+l(«) + « > 0, 

/0,2m+u4j 

(85) /0,2,„(^) + /l,2m(«) + /2,2m(«) = 0, m> 0, 

(86) /0,2m+l(.S) + /l,2m+l(5) + /2,2m+l(5) = 0, m > 0. 

in Laurent series about £, - 1, m > 



(87) /2,2,„(^) = £]^J?l;,,„(.£-l)^, 

CO 

(88) /2,2m+l(.£) = Yj - 1)^' 

J=-l 

(89) fo,2,n{^) = YjCj,,n{^-iy, 

CO 

(90) /o,2m+l 

i=-i 

CO 

(91) /i,2„,(.^) = J]^ £,,„(£ -1)^ 

(92) /i,2,„+i(.£) = 2]^n'«(«-l)'- 

By solving for the leading coefficients of these expansions we can determine the inhomogeneous term 
because 

(93) igL(l) = ei,m^l,m, 

(94) k2,„+l(l) - ^^3,mS_i,™, 
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for m>0. 

Our solution schema is as follows. We first note that coefficient of {£, - of lISTI l gives 

(95) £)_i,„,_i = i^^, m>l. 
Now the coefficient of - 1)° of l|8l]l gives 

(96) — ^ = -£)o,m-i, > 1, 
which combined with (|95] | yields 

(97) Do,,, = -4(m + l):^|^, m>0. 

Next we combine the coefficient of (£ - 1)~^ of ll86l l with the coefficient of {£, - 1)"^ of l l83ll to deduce the 
relations 

4m + 2 

(98) n.i,„, = —- , m>0, 

and 

4m + 2 

(99) S_i,„, = — , m>0. 

Continuing we combine (|95l l and (|99t to derive 

2m + 1 

(100) Ci,,„ = -J?Ii,„,+i, m > 0. 

zm + / 

Now we use the coefficient of (<5 - 1)" of l l83ll to give 

(101) C2,„, = -^^So,,nCl,„, m > 0. 
Proceeding we employ coefficient of {£, - 1)^ of | |83l l to give 

(102) Si,„, = -J?Ii,,„ - Ci,„, + ^^^'"g^^-* [C^ ,„ - Ci,,nC3,,„] , m > 0. 
Proceeding further the coefficient of (£ - 1)^ of lISTIl under m ^ m + 1 yields 

^(fH + 1) r -1 

(103) Di,,„ = -Jli ,m+l — <-'l,m+l + 

Now we employ the result for D-i_m from l|95] l, Di,,;, from 1 11 031 and £)o,m from l|97t to compute that 

,111+1 

(104) £)_i,„,£)i,„, = -4(m + l) 16{m + lY— , m > 0. 

Into the coefficient of (^ - l)'^ of l l84l l we employ the results of l| 10411 and ||95) and compute 

(105) %„, + 2)3,,„ = ^ ^. Jll [^i,,„+i + Ci,„,+i] , m > 0. 

4(m + 1) '-'"'^^ 



On the other hand if we take the coefficient of (cj - 1)'' of I 

(106) S3,,„ + D3,,n + C3,,„ + [-Sl„, + S-l,mSiJ = 0, m > 0, 
and use l IlOll for So,m, l l98l l for S-i,m and (|102l l for Si,,,, we arrive at 

(107) %„,+£>3,», = 2(2;^^[^i,,„+Ci,„], m>0. 
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Combining lll071 l and lIlOSI l, with the assistance of lIlOOI l, gives us a first order difference equation 

Im + 1 



(108) 



l,m+l ,m+l 



2m + 2 

Using the initial data o + Ci,o = -tt"^^^ this solved by 

r(m + i) 



[J?li,,„ + Ci,,„] , m > 0. 



(109) 



m > 0. 



nT{m + 1)' 

Employing l llOOt once more gives us the first order inhomogeneous difference equation 



(110) 



2m + 1 



T(m + i) 



2m + 2" ■^'"'"'^ nT{m + iy 
which is soluble by multiplying through by a summation factor and the end result (J?Ii o = 0) is 



m > 0, 



(111) 



r(m + l) '^' THr+l) 
> , m > 0. 

7ir(m + i)t;^ r2(r + i) 



Back substituting into (|95l98ll00ll , the coefficient of - 1)^ of 
for m > 0, respectively 



and the coefficient of {£, - 1)^ of (|86] l gives 



(112) 

(113) 

(114) 

(115) 
and 
(116) 



2)-l,m - 



4(m + 1) 

^l,m+l 

4(771 + 1) 

^l,m+l 



_ Im + l 
2m + 2 

_£(m+j)_ 
7ir(m + l)' 



1 r(7« + 1) 



47Z (m + l)r(7M + 2) 

Using these formulae and the relations (193194b and |(7T] | gives the evaluation l|80ll . 



By drawing together all the previous results we present the simplest recurrence system for the variance 
of the index. 

Corollary 4. The even and odd sub-cases of the variance of the index satisfy second-order, linear inhomogeneous 
difference equations 



(2m - 1) 



1 + 



and 



ni8^ 

^^^^^ (2^^^"^^ 
respectively. 



1 + 



ilm + 2) (2m + 1) 
(27«)(2m - 1) 

{2m){lm - 1) 



(2m + 2) (2m + 1) 



(2m - 2) 



A 

A2m+1 + t;; ZT^lm-l 



(2m + 2) 



1 r(m + |)r(m-|) 
"27i2 r2(m + l) 

1 r2(m + 1) 
27t2 r(m + 2)r(m + l)' 



m > 1, 



77i > 1. 
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Remark 4. An alternative system of recurrences to that of II20I21I I was discussed in IITTl and was necessary 
in order to relate the system to the first discrete Painleve equation. This system is formulated in terms of 
X2n+i '■= fi.n and X2n+2 := /o,n and the system of recurrences governing these are, for our specialisation 

(119) A:2m(5)[X2m+l(4) + A-2m(.£)+A-2m-l(5)] =2m, WI > 1, 

(120) X2m^l{i)[X2m+2{^)+ X2m+l{C)+ K2m{^)] = '2.ni, M > 0. 

However there does not seem to be any utility in employing this for our problem. 

To conclude this section we give what appears to be the simplest summation representation for the 
variance of the index. 



Proposition 11. Let m = [n/lj. The variance of the index has the summation representations 
n n 



r2(/ + i) 



(121) ^" = i - ^ L ^([TT) ['^('" + 5) - '/'('« + 1) - n + 5) + n + 1)] , 

(122) = ^ |['/'(- + 1) - + i)] g + [pn + 1) - ^[i + i)]| , 

for all n >0. 

Proof. This follows by combining l(77|l and | [80l l, and the functional equation of the digamma fimction \p{x) 
(see 5.5.2 of ||22|| ). The second form is derived using the identity 

-fi T^{1 + 5) r 1 1 n2 

(123) El^[V'(' + i)->/'(' + 5)]-y 



4. Integral Representations 

We now proceed to derive a number of integral representations which while in and of themselves are of 
limited utility they will offer a pathway towards an asymptotic formula. 

Proposition 12. Let m = \ nl2\ with n>0. The variance of the index has the quasi-Beta integral representation 

(124) A„ = ^ -| t-^l\\ - f)-^K( VT^) {2[ip{m + 1) - 4,{m + i)](l - t'") - t'" log(0| , 

where K(f) is the complete elliptic integral of the second kind (see 13.8 o/||9l). 



(125) K(z) = -2fi(l,l;l;z2), 
Proof. We define an auxiliary function 

T^(x + 1) T^{x - m + i + /) 

(126) Rx) := ^ J y -^^ %{x) > -1 m > 0, 

^ ' ' T^x+l)j^T^x-m + l + iy ^'2' 

which serves as a stepping-stone to our goal because of the fact 



n 



1 r2(m + \) d 



^'''^ ^'-4f^^r>7TT)^^ 

Now we can employ the integral representation for the square of the ratio of Gamma fxmctions given by Eq. 
(5.50), p. 198 of (m, which states for fi = a = x- m + ^, 5 = y = x-m + l 

T^(x -m+^ + l) 2 ,n , I 
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This is valid because > + x - m) and ^R(2y - 2a) = 1 > 0. Interchanging the summation and 

integration then gives 

(129) /(x) = l!i^±^i f\te-"'-'/^^f^K(-/T^t), 

r2(x + i)7iJo i-t ' 

where we continue to assume %{x - m) > -j. Employing this into our formula for A„ and the identity 

(130) dt ri/2(i - iog(f )K( Vr^) = - y , 

we deduce H124b . □ 
Let us define a function similar to that of l|126|l 

(131) /,„:= dtt-'/^——K{^^n), %{m) > 0. 

Jo ^ ~ ' 

One can then verify 

(132) A„ = ^[^{m + 1) - + l)]/„, + ^^/m 
holds where m - [n/l] and n > by using the definition of /„, and (I124I I. 

Proposition 13. Again let m = [_nl2\ and n> 0. The variance of the index has either of the Mellin-Barnes integral 
representations 



(133) A„ = — J^_^^ — 2 [^(m + 1) - ^(m + 1)] ^^^^ [^(m + s) - ^(s)] 



r2(m + 1 - s) 

ds r2(i - s) 



or 

(134) 

^ r2(l_s) + 1) - + !)]['/'('« + s) - •/'(s)] + [i/'(i - s) - V^(l - s)]!/'(m + s)} , 

where < c < |. 

Proof. From the tables of fT9ll we find two inverse Mellin transform formulae 

(135) — — = —X-' [ii^{m + s) - ip{s)] , c> 0, ^}l(w2) > 0, 

X — L Jc—icci ^Tll 

log(x) r^"" ds 

(136) c>0. 

We insert these into l|1241 and interchange the orders of the integration and are left with the following type 
of integral, which is evaluated using Eq. (14.2), p. 155 of ||T91 

(137) r rftf-^K(Vr^) = ^ ^^^^\ , ^)l(s)>0. 

Jo 2 r2(s + i) 

By employing this integral we are obliged to keep c within the analyticity strips (0, |) and (0, m-V }f) for the 
first and second terms respectively. This yields ( |133l l whereas to deduce l |134t we integrate the second term 
of (|133|l by parts (which is actually valid for 0<c<m + i) and displace the contour under s i-^ s + 77i. □ 

Remark 5. Employing the standard Mellin transform techniques as expounded in, say [231, to extract the 
n — > oo asymptotics from either i ll 331 or 1 11 341 does not yield a clean large n expansion. 

We give a third integral representation which bears some resemblance to the Schwinger proper- time inte- 
gral representations in Quantum Electrodynamics of particle-antiparticle systems subject to large imiform 
backgroimd electromagnetic fields, see Eq. (5.25) of @ or Eq. (7.88) of (121. 
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Proposition 14. The variance of the index has an integral representation on (0, oo) with 

Jr^ca rl { (-\ 1 

— ^ J 1 _ cosh-2'" X2Fi(i + 1; cosh-2 x)\ , m>0. 
sinhx m! ^ 

Proof. We employ some elementary changes of variable and the elliptic integral definition of K(z) in order 
to rewrite the integral as 

dy^= dx {x^ + y^r^ l-(l + x^ + 
CO yl + y"^ J-oo L 

We evaluate the inner integral in terms of the Gauss hypergeometric function 'Kiy) i= 0,m e Z, m > 

and employ the connection formula Eq. (15.10.29) of 1122 1 to rewrite the hypergeometric function. The final 
result is ( IBSl . □ 

5. ASYMPTOTICS AS n ^ CO 

None of the above integral representations yield the asymptotic behaviour of the variance of the index 
as « ^ 00 in a particularly clean way and we therefore present an evaluation which while being exact also 
constitutes an asymptotic formula. 

Proposition 15. The function }„, has the evaluation 

(141) /,„ = Itt [^{m + i) - ^.(1)] + in log(4) - ^^ ^(1. h I h + 1)' 
for 'K{m) > 0. Then the variance of the index is given by ( I132P . 

Proof. Into the definition (|131l l we employ the series expansion for K(z) about z = 0, interchange the 
integration and summation and integrate term-by-term after separating out the first term from the rest. For 
the first term we employ the direct Mellin transform corresponding to (|135|l and the standard beta-integral 
for the remainder. In this result we recognise the series expansions of the hypergeometric functions 3F2 and 
4F3 with argument unity and recall the identity 

(142) 3f2(i,l,l;2,2;l) = 21og(4). 

This gives the expression 1 11411 . □ 

Clearly lll32b and (I141l l give the large n behaviour in an explicit way with the leading logarithmic term 
coming from !/'(m + |), which we detail in the following corollary by giving the leading order terms. 

Corollary 5. The leading order terms of the asymptotic expansions ofA„ asn 00 are, in the even case n = 2k 



log 8k + y + 1 log 8k + y 7 
(143) A21, = ^ ^ / ° ' + 



+ 



27t2 8n^k 38An^B 

241og8fc + 247-41 219 2560 log 8fc + 2560y - 6247 19129 



1536Ti2fc3 819207x2^4 3276807i2fc5 16515072n^k^ 

and the odd case n = 2k+ 1 

logSfc + y + l log8/c + y + 2 24 log 8/c + 24y - 7 
(144) A2,.i = + 

8 log 8A: + 8y - 9 1920 log 8k + 1920]/ - 3937 
5127i2fc3 2457607i2fc4 
2560 log 8k + 2560y - 5809 322560 log 8k + 322560y - 882767 
3276807T2fc5 825753607i2fc6 
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+ o(r7iogfc). 



■ o(r^ log fc). 



The asymptotics for the variance of the index can also be studied via the Birkhoff-Trjitzinsky theory 0, 
Q/ ID/ lEI for linear difference equations. However this is unnecessary in our case because of the following 
fact. The homogeneous forms of the linear second-order difference equations Illl7lll8t are 



and 



1 + 



(2m + 2)(2m + 1) 



1 + 



(2m)(2;7i - 1) 
(2m)(2OT - 1) 



{2m + 2){2m + 1) 



(2m) _ 

+ (2^;r:2)^'"-^ ~ °' 



(2m) _ 
^'"+(2^;^^'"-^-°' 



which have the general solutions 

(147) X„, = Ci(m);7z + C2(m)m [log 4 + i/;(m + i) - i/'(m + 1)] 

(148) X„, = Ci(m)(m + i) + C2(m)(m + i) [log 4 + i/;(m + i) - i/;(m + 1)] , 

respectively, for arbitrary periodic functions C\{m), C2(m) with period unity. Thus the asymptotic behaviour 
is clearly apparent in our simple solutions. 
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